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Abstract 

We derive an analytical expression of the inter-arrival time distri- 
bution for a non- homogeneous Poisson process (NHPP). This expres- 
sion is exact and is applicable to any time interval, finite or infinite. 
As an illustration, we present simulation results for three different 
intensity functions. 

The Poisson process has found numerous applications in science, engi- 
neering, economics and other areas. The NHPP is probably the best known 
generalization of the Poisson process (see for example [1]). It is character- 
ized by a deterministic intensity function that describes how the rate of the 
process changes in time. For an ordinary Poisson processes, this function 
is a constant. The NHPP has many important characteristics, but in this 
work we concentrate on the distribution of inter-arrival times, also known as 
a waiting time distribution, a somewhat "cumulative" property of the pro- 
cess. For an ordinary Poisson process, this distribution is an exponential if 
the process is observed over a sufficiently large time interval. In the general 
case of the NHPP, various approximations were developed; examples include 
the time-dependent Poisson process, when intensity changes gradually with 
time, and the conditional Poisson process, when the intensity function has a 
probabilistic interpretation. 

Let us state the central question of this work: Given that the NHPP starts 
at time and that we observe events during a fixed (finite or infinite) inter- 
val T, what is the distribution of inter-arrival times between consecutively 
arriving events? 
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We denote the arrival times of the NHPP by to = 0, ti, . . . , t„. The 
corresponding inter-arrival times are Si = ti, S2 — t2—ti, . . . , Sn — tn—tn-i- 
Let \(t) (intensity function) be a known, deterministic function of time and 
let A{t) = £\{u)du < 00 be the mean value function (expected number of 
events of the NHPP on the time interval (0,t)). The total expected number 
of events in the interval T is A(T). Let N{t) be the number of arrivals on the 
interval (0, t). The probability of having n arrivals in the interval {t, t + x) 
is given by [1] 

Prob{N{t + x)- N{t) = n} = [A(t + ^) - A(^)]\ -[A(t+.)-A(t)] 

nl 

where conventionally we set 0! = 1. 

First we write an equation for the probability of the first inter-arrival 
time si > X and the total number of events N{T) — n 

A{x,n) = Prob{si > a; p| N{T) = n} 

_ [A(r) - A(a;)]" e'^^^) 
~ n\ A(T) 

where 1/A(T) is a normalization constant (inverse of the expected number of 
events on the interval T), exp(— A(a;)) is the probability to have no arrivals 
on the interval (0, x) and the second factor is the probability to have exactly 
n arrivals on the interval (x, T). These probabilities should be multiphed, as 
N{b — a) and N{d — c) are independent random variables on disjoint intervals 
(a, b) and (c, d) [1]. 

Next, wc write an expression for the probability of the A;*^ arrival time 
tk E {y,y + dy), k > and the A;*'* -I- 1 inter-arrival time Sk+i > x and the 
total number of events N{T) — n 

dB{x, y, k, n) = Prob{tk E iy,y + dy) Q Sk+i > ^ N{T) = n} 
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where the first factor is the probability to have the k arrival time tfc E 
{y, y + dy) [1], the second factor is the probability to have no arrivals on the 
interval {y,y + x), and the last is the probability to have n — k arrivals on 
the interval {y + x,T). Again, these probabilities should be multiplied due 
to mutual independence. 

The total probability to have an inter-arrival time S > x, G{x), is given 

by 

oo oo n— 1 pT—x 

G{x) = Prob{S >x} = J2 ^(^' ^) + ^Yl / ^' ^' ^) 

n=l n=2 k=l ^ 

where the upper limit of integration is T — a;, because for any given x < T, 
y + x must be less then T. Interchanging the order of integration and the two 
summations {X{t) is integrable function), summing up by k and n, integrating 
over y and by parts, we arrive at 

T-x 

A(y + x)e-[^(^+")-^(^)ldy (5) 

We see that this probability is properly normalized, G{0) = 1 and G{T) = 0. 
The probability density function (pdf) for Equation (0) is given by —G^{x) 
and (after integrating by parts) is equal to 

pdfix) = ^ \{y)\iy + x)e-[^(^+-)-^(^)]rfi/ + ^^^^ (6) 

This equation can be interpreted as follows: The factor X{y + x)exp{ — [A{y + 
x) — A(?/)]} is the pdf to have the inter-arrival time x given that the sum 
of all previous Sj is y. \{y) / K{T)dt is the probability for the NHPP to be 
found in the interval (y, y + dy), y > 0. This means that at least one arrival 
of NHPP has already happened. In other words, the first arrival density is 
missing in the integral, which is exactly compensated by the second term in 
Equation (0). 

To illustrate Equation © we perform numerical simulations for three 
intensity functions A(t): X{t) = 1/c for an ordinary Poisson process, A(t) = 
c/(l -|- t) for a decaying intensity function, and A(t) = ct for a growing 
intensity function. These functions were chosen because the integrals in 

can be taken explicitly. Our simulations were done as follows: For every 
intensity function we ran 10^ processes. During every process, the next inter- 
arrival time Sj+i was drawn from a known function (via the Inverse Transform 
method for NHPP, see [1]) which depends on the sum of all previous Sj's. 
This inter-arrival time was utilized for statistics unless the corresponding 
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Figure 1: pdf for X(t) = 1/c. Solid lines are from Equation ((Tj) and dashed 
lines are for the corresponding exponential distributions. 

arrival time tj+i > T, in which case the current process was terminated and 
the next was started. 

On all our figures simulation results are presented in colors while the- 
oretical pdf's from ^ are given by solid black lines. The distribution of 
inter-arrival times for X{t) = 1/c is given by 



and examples are shown in Figure (1) for three different values of c. This is a 
homogeneous Poisson process with the mean time between events c. Dashed 
lines represent three corresponding exponential distributions. The deviations 
from an exponential distribution are easily understood as follows: given that 
the process is observed at time y < T (arrival time ti = y), the inter-arrivals 
time Sj+i > T — y cannot be observed because the corresponding arrival times 
will exceed the observational interval T. In other words, the finiteness of T 
introduces memory into the process. 

Our next example is for the intensity function X{t) = c/(l + t). The 
corresponding inter-arrival pdf is given by 




(7) 
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Figure 2: pdf for X{t) = 1/c. Solid lines are from Equation (jSJ. 

and simulation results are shown in Figure (2). 

The last intensity function to be considered here is X{t) = ct. Simulation 
results are presented in Figure (3). An explicit expression for the inter-arrival 
pdf can be easily found using Mathematica and is not presented here because 
of it's length. 

In this work we obtained an exact expression for the inter-arrival time 
distribution for NHPP. Results of this work will be applied to the distribution 
of earthquake aftershock waiting times [2]. 
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